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Abstract 

We  propose  a  new  approach  for  designing  diffractive  optical  corrective  elements  with 
zooming  capability  to  convert  nonlinear  sinusoidal  scanning  into  linear  scanning.  Such 
devices  will  be  useful  for  linearizing  the  angular  scan  of  a  resonant  mirror  scanner.  The 
design  methodology  is  to  create  a  graded  index  of  refraction  device  as  the  reference 
design,  with  its  index  of  refraction  parameters  based  on  beam  retardation  through 
propagation  in  an  inhomogeneous  media.  In  a  new  approach,  we  make  explicit  use  of  the 
symmetry  of  the  optical  problem  to  mathematically  derive  the  form  of  the  refractive 
gradient.  The  diffractive  element  is  designed  utilizing  a  binarizing  algorithm  of  the 
accumulated  phase  from  transmission  through  the  refractive  element. 


1.  Introduction 


Diffractive  optical  elements  have  found  considerable  use  in  many  areas  of  optics.  One 
important  area  of  application  is  as  a  corrector  element  in  optical  systems  to  modify 
wavefronts  to  remove  distortions  or  aberrations.  For  example,  a  diffractive  element  has 
been  designed  to  provide  compensation  for  thermally  induced  aberrations  in  a  solid-state 
laser1.  In  applications  such  as  these,  diffractive  technology  provides  a  planar, 
lithographically  fabricable  solution  to  replace  difficult  to  realize  3-D  conventional 
refractive  optics  .  Diffractive  elements  can  also  be  utilized  in  important  applications  such 
as  beam  steering3  and  optical  free  space  interconnects4.  In  the  present  work  we  seek  to 
find  a  diffractive  solution  to  the  important  practical  problem  of  distortion  in  fast  optical 
scanners. 

High  speed  optical  scanning  of  laser  beams  has  found  numerous  applications  in  science, 
engineering  and  medicine,  from  the  ubiquitous  laser  printer  to  geophysical  remote 
sensing  and  confocal  microscopy.  An  optomechanical  device  known  as  a  resonant 
scanner  is  often  used  for  scanning  lasers.  Resonant  scanners,  which  are  oscillating 
mirrors  excited  to  operate  at  their  natural  resonant  frequency,  provide  a  nonlinear 
sinusoidal  scan,  and  at  the  limits  of  the  scan,  where  the  sine  function  changes  more 
slowly  near  its  peak,  the  scan  field  is  significantly  compressed.  Despite  this  deficiency, 
resonant  scanners  have  many  attractive  attributes  and  are  widely  used.  In  imaging 
applications,  the  scanned  image  can  be  electronically  corrected  during  or  after  acquisition 
to  produce  an  image  equivalent  to  which  would  be  produced  by  an  ideal  linear  scanner,5’6 
whereas  other  applications  require  that  the  scan  itself  must  be  made  linear.  An  optical 
means  of  correction  would  be  highly  desirable,  since  it  means  that  the  wide  variety  of 
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commercially  available  resonant  scanners  could  be  utilized  without  any  modification  to 
the  scanner  itself. 

In  previous  work,  we  have  reported  on  the  derivation  of  a  diffractive  element  and  its 
associated  refractive  element  for  the  very  same  purpose.  The  derivation  of  the  graded 
index  of  refraction  of  the  refractive  active  element  in  our  prior  work  was  based  on  beam 
propagation  in  an  inhomogeneous  media  and  on  the  beam  retardation  approach  .  The 
diffractive  element  is  simply  a  binary  phase  or  amplitude  version  of  the  refractive 
element.  In  the  present  work,  we  describe  a  new  design  based  on  the  beam  retardation 
approach  (phase  lag),  which  explicitly  takes  into  account  the  symmetry  of  the  optical 
problem.  As  with  our  earlier  work9,  which  we  will  refer  to  as  the  ‘asymmetric  solution’ 
the  present  design  does  take  the  beam  diameter  into  consideration. 

Comparison  of  the  first  order  approximation  solutions  of  the  beam  propagation  and  the 
beam  retardation  approaches,  suggests  that  for  a  quasi  broad  beam,  the  solution  from  the 
beam  propagation  approach  should  collapse  to  the  solution  from  the  beam  retardation 
approach  at  this  level  of  approximation.  As  we  shall  see  later,  the  differential  equations 
we  develop  only  contain  a  linear  dependence  on  the  ratio  of  the  element  size  to  the  beam 
diameter  AD.  These  considerations  suggest  that  this  new  variant  of  the  beam  retardation 
approach  will  work  well  with  narrow  to  moderately  broad  beams.  The  addition  of  the 
higher  order  terms  will  naturally  improve  the  accuracy  of  our  computation. 

2.  Device  Theory 

Our  theory  for  developing  the  index  of  refraction  of  the  corrective  element  for  converting 
sinusoidal  to  linear  scanning  consist  as  follow:  In  Fig.  1,  let  us  assume  that  we  have  a 
corrective  optical  element  which  is  located  at  a  distance  equal  to  Z2  from  a  sinusoidal 
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scanner.  The  reflected  beam  out  of  the  scanner  is  incident  upon  the  corrective  element  in 
the  y  axis  at  I12  and  at  an  angle  02  which  is  proportional  to  sin  (cot),  where  co  is  the 
scanning  angular  frequency.  The  corrective  element  deflects  the  output  beam,  which 
propagates  at  angle  63,  and  meets  the  target  plane  at  point  b.  This  point  b  is  located  at  the 
target  plane  at  the  highest  of  h  from  the  axis  and  at  a  distance  of  Zi-Z2  from  the  corrective 
element  (or  Zi  from  the  scanner).  Point  b  corresponds  to  the  point  where  the  beam  should 
arrive  if  the  sinusoidal  scanner  is  replaced  by  a  linear  scanner.  The  corresponding  virtual 
scan  angle  is  9 1,  as  shown  in  the  diagram. 

For  developing  the  design  of  the  required  corrective  element,  we  use  the  beam  retardation 
approach.  In  this  method,  we  consider  that  we  have  a  corrective  element  with  a  graded 
index  of  refraction.  When  a  beam  propagates  through  the  graded  index  of  refraction 
material,  it  causes  one  side  of  the  beam  to  phase  lag  behind  the  other.  This  phase  lag 
causes  the  beam  wavefront  to  bend.  We  derived  the  exact  formula  for  this  refractive 
element;  however,  since  it  is  hard  to  fabricate  a  refractive  element  with  that  level  of 
complexity  of  the  gradient  in  the  index  of  refraction,  or  its  corresponding  diffractive 
element,  we  use  a  binarization  scheme  to  modify  the  results  to  produce  a  diffractive 
element  that  is  easier  to  fabricate. 

Before  we  derive  the  fonnula  for  the  corrective  element,  it  is  essential  that  we  derive  the 
output  angle  #3  from  this  optical  element  as  a  function  of  B\,  the  desired  linear  scan  angle. 
That  is,  we  want  a  closed  form  mathematical  expression  of  the  functionality  of  the 
element. 

From  simple  geometrical  considerations  in  Fig.  1  it  is  possible  to  show  that  the  linear 
scan  angle  9i  and  the  input  scan  angle  02  to  the  corrective  element  are  given  by, 
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(1) 


h  =  Zj  tan  Ox , 
and, 


h2  =  Z2  tan  02 .  (2) 

From  Fig.  1  and  equations  1  and  2  it  is  possible  to  show  that  the  output  deflection  angle 
03  is  given  by, 


Tan  0,  = 


7  -7 

A  /-j  2 


7,  tan  0]  -  Z2  tan  02 
7  -7 

A  /-j  2 


And  in  the  paraxial  approximation, 


(3A) 


03  = 


h-h2  zxex  -z2e2 


7-7 

z^i  Z,2 


7  -7 

Z.]  Zv2 


(3B) 


For  sinusoidal  scanning  it  possible  to  express  the  scan  angle  02  in  the  form, 


y 

02  =  02O  sin  cot  =  — , 


(4A) 


where  6bo  is  the  maximum  scanning  amplitude  and  co  is  angular  frequency  for  scanning. 
For  calculating  the  graded  index  of  refraction  of  the  refractive  reference  element,  as  a 
function  of  the  y  coordinate,  it  is  necessary  to  express  the  input  and  the  output  angles  in 
tenns  of  the  y  coordinate  instead  of  the  time.  Equation  (4A)  can  be  written  in  the  inverse 
function  form  as: 


t  =  —  suT1 — — —  (4B) 

co  Z202{) 

On  the  other  hand,  if  we  assume  that  the  sinusoidal  scanner  is  replaced  by  a  linear 
scanner,  then  it  is  possible  to  write  the  equation  of  motion  as 

0i  =  gO  (5) 

where  g  is  a  constant.  Substituting  the  time  value  of  equation  4B  in  equation  5  yields 
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0,  =  —  sin  1  -y- 


CO  l  ^2^20 


Substituting  equations  6  and  4 A  in  equation  3,  yields  the  output  of  the  diffracted  beam  as 
a  function  of  y,  given  by  the  following  equation: 


Z,  —  sin  1  -OL 


tan  8,  «  8,  = 


-z,i-  gzx sin-1  -y- 
Z2820)  Z2  { Z2820 


(Zi-Z2) 


(Z\  Z2  )c 


which  mathematically  represents  the  action  of  our  corrective  element.  We  now  proceed 
to  the  detailed  design  of  the  element. 

2.1  Beam  retardation  approach 

Let  us  assume  as  shown  in  figure  2  that  there  is  a  beam  of  diameter  AD  propagating  in  a 
graded  index  of  refraction  medium  of  thickness  /.  The  medium  has  a  graded  index  of 
refraction  in  the  y  direction.  As  mentioned  earlier,  during  propagation,  one  side  of  the 
beam  phase  lags  behind  the  other  side,  which  in  turn  makes  the  beam  wave  front  bend. 
The  phase  lags  of  the  upper  and  lower  sides  of  the  beam  respectively  are: 

„  ln(8h  +  A»)/ 


and  d5  = 


2x(n0)l 


Therefore  the  beam  wavefront  bends  by  an  angle  of 


A  645  =  #4  -  #5  = 


2;r(A  n)l 


The  deviation  as  a  result  of  the  phase  lag  can  be  written  as, 


A#45  = — An(AD)£ . 
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For  conversion  of  sinusoidal  scanning  to  linear  scanning,  the  beam  deflection  due  to 
propagation  through  the  refractive  element  should  be  equal  to  the  difference  in  the 
propagation  angle  between  sinusoidal  and  linear  scanning,  e.g. 


A#45  -  A $32  -  $3  $2  • 


(12) 


Equalizing  equations  1 1  and  12,  and  using  equation  3B  for  the  value  of  6%  yields 

(13) 


—  An(AD)£  =  $3  -e2  =  Z'°'  Z2°2  -  0, . 


A 


7  -7 

Z>1  2 


Substituting  the  values  of  6\,  $2  from  equations  4A  and  6  into  equation  (13)  yields 

(14) 


— An(AD)l  =  ^sin _1 
A  co 


'_y_' 

VZ2^20  J 


-A*-. 


Since  a  resonant  scanner  deflects  symmetrically  about  its  central  rest  position,  the 
corrective  optical  element  that  corrects  the  deflected  beam  for  nonlinear  scanner 
deflection  must  have  a  graded  refractive  index  that  is  symmetrical  with  respect  to  the 
central  or  zero  position. 

In  order  to  solve  the  differential  equation  (14),  we  need  to  replace  An  by  its  Taylor  series 
expansion.  For  a  symmetric  function  around  the  origin,  An  can  be  reduced  further.  The 
Taylor  series  expansion  of  the  index  of  refraction  at  the  origin  is, 


/  \  dn  d2n  2  <3 3n  3  d4n  4 

i[x  }  =  nf]+  a, — x  +  a7 — r-x  +0, — rx  +a. — rx  +■ 
ox  ox  ox  ox 


(15) 


where  n0  is  the  offset  (the  refractive  index  at  the  center),  and  ai,  a2,  a3  etc.  are  given 

by  a,  =  1  ,  a,  =  —  ,a,  =  —  and  so  forth. 

2!  3!  4! 
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Since  the  diffractive  element  index  of  refraction  n(x)  should  be  symmetric  around  the 
point  of  symmetry,  the  terms  containing  the  odd  powers  in  x  should  be  eliminated, 
thereby  reducing  the  Taylor  series  expression  for  the  index  of  refraction  to: 


/  \  d2n  2  d4n  4 

n{x)  =  n0  +  a2—Yx  +a4—Tx  +•••. 
ox  ox 


(16) 


The  graded  index  of  refraction  An(x)  can  now  be  derived,  and  therefore, 

^2  ^4  ^  6 

A  ,  .  „  o  n  .  .  3  o  n  K  r  5  c  n  K 

An(x)  =  2a2x — -Ax  +  4a4x  ^-j-Ax  +  6a6x  ^-^-AxH - . 

dx 


8xq 


dxb 


(17) 


Substituting  expression  (17)  for  An(x)  in  equation  (14),  and  remembering  z2&io=yo, 
yields, 


2  n 

T 


r\2  ^4  *~\6 

on  A  q  o  n  ^  5  on 


2a2y—^  +  4a4y  — 
oy  ay 


+  6  a6y 


dyb 


(AD)/  =  ^sin  1 


<x> 


vTo  j 


-A*-. 


(18) 


where  a 2  =  — ,a4  =—  ,a6  =-. 
2  2!  4  4!  6  6! 


It  was  not  possible  to  directly  solve  differential  equation  (18).  To  overcome  the 
problem  of  not  having  a  direct  solution,  the  arcsine  on  the  right  hand  side  of  the 
equation  was  replaced  by  its  polynomial  expansion  to  the  seventh  order.  Figure  3 
shows  sin  (x)  and  its  polynomial  expansion  to  the  seventh  and  the  ninth  orders.  Using 
this  approach,  Eqn.  (18)  was  solved  for  both  the  second  and  the  fourth  order  Taylor 
expansions  of  the  index  of  refraction. 

The  solution  for  the  second  order  expansion  is, 


n{y)  = 


5  B0y0 


6272 


_y_ 

y0J 


+ 


At0 


400 


z_ 

To 


+ 


At0 


72 


_r_ 

To 


+ 


(A)  -£|To)tq 


_T_ 
To  J 


+  C],  (19) 
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where  B0  = 


AgA 

2n(tyD)co\ 


,  and  B]  = 


AA 

2x(A D)Z2l ' 


Similarly  the  solution  for  the  fourth  order  expansion  is, 

5B0y0 


n(y ) : 


225792 


r  \ 

y 

s  +  B«y0 

f  \ 

y 

6  +  B0v0 

f  \ 

y 

+  y2(Bo 

~6Bi  y0)y0 1 

y 

UJ 

7200 

UoJ 

576 

UoJ 

12 

UoJ 

+  Q- 


(20) 


The  details  on  the  derivation  of  the  above  equation  can  be  found  in  the  appendix  A. 

The  ratio  g/to  in  both  equations  14  and  18,  which  is  also  expressed  within  the  parameter 
Bo  of  equations  19  and  20,  has  significance  in  the  design  of  the  diffractive  element.  The 
ratio  of  g/to  determines  the  zooming  capability  or  the  magnification  of  the  corrective 
element  either  to  be  equal,  less  or  larger  than  the  original  uncorrected  field  of  view. 
Looking  at  Fig.  4,  which  compares  the  functions  of  motion  for  both  linear  and  nonlinear 
scanning,  reveals  the  following:  in  order  for  the  corrective  element  to  have  no  effect  on 
the  dimensions  of  the  field  of  view  (FOV),  it  is  necessary  for  the  amplitudes  of  both  the 
sinusoidal  and  triangular  motions  to  be  equal.  This  is  satisfied  when  g  =  2(o/n  and  g/to  = 
2/71.  On  the  other  hand,  in  order  for  the  diffractive  element  to  maintain  the  linear  portion 
of  the  sinusoidal  scanning  near  the  nodal  line,  the  tangent  of  the  triangular  function 
should  be  equal  to  tangent  of  the  sinusoidal  function;  e.g.  g  =  co.  Also,  an  examination  of 
equations  19  and  20  shows  that  increasing  the  ratio  yo/AD  (remembering  that  B0  and  Bi 
contain  1/AD)  will  require  an  increase  in  the  index  of  refraction  n  for  the  corrective 
element. 

2.2  Binarization  Algorithm 

Fabrication  of  a  corrective  element  with  a  complex  graded  index  profile  based  on 
equations  19  and  20  is  not  easy.  The  gradient  index  refractive  element  is  therefore  best 
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used  as  a  reference  design.  Hence,  we  propose  to  build  a  binary  diffractive  element  for 
image  corrections  that  has  the  following  transmissivity: 


For  binary  amplitude  encoding: 


T  =  \ 


An  A  TC  ATE 

2 (m  - 1  )—<—«/<  2m  —  =  0 
b,  A  b, 

An  An  An 

(2m  - 1)  —  <  —  nl  <  (2m  + 1)  — 

b,  A  b, 


=  1 


(21) 


and  for  binary  phase  encoding: 


T  =  \ 


An  An  An 

2(m-X)——<——nl  <2m——  =  -\ 
b,  A  b, 

(2m -\)—<—nl<  (2m  + 1)  —  =  1 
b,  A  b, 


(22) 


where  b[  is  the  number  of  binarization  levels;  if  bi  =2,  the  encoding  is  binary,  bi  =3,  the 
encoding  is  ternary,  if  bi  =4,  the  encoding  is  quaternary  and  so  on. 

Selection  of  the  offset  value  [the  constant  Cl  in  equations  (19)  and  (20)]  is  extremely 
important  for  the  design  of  the  binary  amplitude  version  of  the  diffractive  element.  This 
is  not  so  critical  in  the  binary  phase  version  of  the  diffractive  element  and  for  the 
refractive  element.  In  the  design  of  the  binary  amplitude  version  of  the  diffractive 
element  one  should  be  sure  that  the  center  of  the  diffractive  element  is  transmittive. 

In  order  to  have  a  transmissive  diffractive  element  at  the  origin,  the  phase  delay  through 
the  transmission  at  the  origin  of  the  refractive  element  should  satisfy  the  following 
condition  for  its  value: 

2m(0)l  /  \2  n 

- ±J-  =  {2m  +  l)— ,  (23) 

A  b, 

where  n(0)  is  the  index  of  refraction  at  the  origin  (e.g.  at  y=0),  and  in  equations  (19)  and 
(20)  is  equal  to  the  constant  of  integration  “Ci”.  Rearranging  (23)  we  have, 
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Ill  =  — 
2 


n(0)lb, 

A 


-1 


(24) 


Setting  n  (0)»1,  is  a  reasonable  selection  of  the  graded  index  of  refraction  at  the  origin,  in 
order  to  estimate  m.  We  thus  obtain, 


1 

m  =  — 
2 


lbL_  i 
/l 


For  the  generation  of  the  diffractive  element  we  selected: 
m  =  round 


A 


v 


and  the  constant  of  integration  from  equations  19  and  20  is  thus, 


2  round 

“li 

fAb 

+  1 

/l 

_2 

l  ^  V. 

lb, 


(25) 


(26) 


(27) 


2.3  Simulations 

We  utilize  commercial  software  (MATLAB,  MathWorks,  Natick,  MA)  to  compute, 
process,  and  display  the  results  obtained  from  our  implementation  of  our  corrective 
element  design.  Figure  5  depicts  drawings  showing  the  diffractive  elements  generated 
from  the  second  order  approximation  of  our  algorithm  by  binarization,  for  a  constant 
aperture  of  1  cm  for  different  beam  diameters  (a)  AD  =  0.01  cm,  (b)  AD  =  0.005  cm  and 
(c)  AD  =  0.002  cm.  From  these  drawings,  we  can  see  that  the  binary  diffractive  element 
consists  of  a  chirped  grating  whose  frequency  increases  as  it  goes  to  the  lateral  edges  of 
the  diffractive  element.  An  increase  in  the  grating  frequency  is  associated  with  an 
increase  in  the  graded  index  of  refraction  of  the  corresponding  refractive  element.  As 
expected  intuitively,  smaller  beam  diameters  require  smaller  grating  spacing  for  the 


11 


diffractive  element  to  function  properly,  and  a  more  exact  reason  may  be  found  by 
examining  Figs.  6  and  8.  Figure  6  shows  plots  of  the  graded  index  of  refraction  for  the 
second  order  approximation  using  equation  17.  In  all  of  these  plots  it  was  assumed  that 
the  refractive  element  thickness  1  =lmm.  (a)  The  diffractive  element  aperture  size  was 
kept  constant  (yQ=  1  cm),  while  three  different  values  of  the  beam  diameter  AD  were 
considered;  (b)  the  beam  diameter  AD  of  the  diffractive  element  was  kept  constant  at  0.01 
cm  while  three  different  values  of  the  diffractive  element  aperture  y0,  0.5,  1  and  2cm, 
were  considered.  It  is  evident  from  the  plots  in  6(a),  that  beams  with  smaller  diameters 
require  a  diffractive  element  with  a  higher  graded  index  of  refraction  in  order  to  achieve 
the  same  beam  propagation  correction,  which  is  reflected  in  the  higher  grating 
frequencies  we  observe  in  Fig.  5.  The  plots  in  6(b)  reveal  that  larger  corrective  element 
apertures  also  require  a  larger  graded  index  of  refraction, 

Figure  7  shows  the  same  plots  in  the  same  order  as  in  figure  5,  but  using  the  graded  index 
of  refraction  for  fourth  order  approximations  derived  in  equation  20.  The  drawings  are  for 
a  constant  aperture  of  1  cm,  for  different  beam  diameters;  (a)  AD  =  0.001  cm,  (b)  AD  = 
0.0005  cm  and  (c)  AD  =  0.0002  cm.  Similarly  Fig.  8  shows  plots  similar  to  those  in  Fig. 
6,  but  derived  from  the  more  accurate  fourth  order  approximations  instead  of  the  second 
order  expressions  utilized  for  Fig.  6.  For  Fig.  8(a),  the  aperture  size  was  kept  constant, 
while  the  beam  diameter  took  on  three  different  values  of  0.02  mm  (0.002  cm),  0.05  mm 
(0.005  cm)  and  0.1  mm  (0.01  cm).  Again,  we  observe  that  beams  with  smaller  diameters 
are  associated  with  higher  graded  indices  of  refraction.  This  can  be  explained  by 
considering  that  a  beam  with  a  larger  diameter  requires  less  graded  index  of  refraction  in 
order  to  achieve  the  same  phase  lag  between  the  two  sides  of  the  beam.  For  a  beam  with 
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an  infinitesimal  diameter,  the  graded  index  of  refraction  needs  to  be  extremely  large  to 
accumulate  a  sufficient  beam  lag  between  the  two  sides  of  the  beam.  In  Fig.  8(b),  two 
different  values  of  the  aperture  size  y0,  1  and  0.5  cm  were  utilized,  whereas  the  beam 
diameter  was  kept  constant  at  0.01  mm  (0.001  cm).  The  larger  aperture  is  associated  with 
a  small  increase  in  the  graded  index  of  refraction. 

2.4  Results  and  Discussion 

We  have  presented  a  formula  for  the  graded  index  of  refraction,  in  equation  17,  and  two 
solutions  of  different  degrees  of  approximation  in  eqns.  19  and  20.  In  each  case,  the  value 
of  the  refractive  index  exhibits  a  linear  dependence  on  the  ratio  of  the  dimension  of  the 
diffractive  element  and  the  value  of  the  beam  diameter  AD,  i.e.  the  use  of  a  refractive 
element  or  its  associated  diffractive  element  fabricated  in  accordance  with  the  appropriate 
value  of  the  beam  diameter  is  thus  desirable,  while  using  a  diffractive  element  intended 
for  correcting  a  different  diameter  beam  is  less  deleterious  than  in  most  other  design 
methods. 

When  we  examine  a  section  of  the  diffractive  elements  shown  in  Fig.  5  or  Fig.  7  we  see 
that  any  portion  of  the  diffractive  element  consists  of  a  chirped  grating.  Thus  we  expect  a 
problem  with  beam  fanning.  If  the  grating  consists  of  equally  spaced  lines,  then  the 
diffraction  pattern  from  such  a  grating  will  be  distinguished  by  diffraction  order  numbers. 
However,  because  this  grating  has  no  constant  spacing,  the  left  side  of  beam  will  diffract 
more  than  the  right  hand  side  of  the  beam  and  this  causes  beam  fanning  in  addition  to 
producing  regular  diffraction  orders  around  the  average  line  spacing.  This  can  be 
minimized  by  careful  attention  to  system  requirements  and  designing  for  the  appropriate 
beam  diameter. 
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3.  Conclusions 


In  conclusion,  we  have  proposed  a  refractive  corrective  element  with  zooming  capability 
for  converting  nonlinear  sinusoidal  scanning  into  linear  scanning.  We  derived 
expressions  for  the  graded  index  of  refraction  for  the  refractive  element  based  on  the 
beam  retardation  approach  through  propagation  in  inhomogeneous  media.  In  a  new 
approach,  we  make  explicit  use  of  the  symmetry  of  the  optical  problem  to  mathematically 
derive  the  form  of  the  refractive  gradient.  An  algorithm  for  converting  the  graded  index 
of  refraction  design  to  the  binary  diffractive  version  is  also  presented.  In  our  theory,  we 
found  that  the  index  of  refraction  increases  as  we  moved  to  the  lateral  edges  of  the 
corrective  element,  behaving  like  a  diverging  lens.  In  the  binary  version  of  this  element, 
this  has  been  expressed  by  having  a  chirped  grating  with  an  increasing  frequency  as  one 
moves  from  the  center  to  the  sides.  We  also  find  from  the  theoretical  derivation,  and  as 
expressed  in  the  computer  simulations,  that  the  graded  index  of  refraction  increases  with 
the  requirement  of  an  increase  in  the  aperture  of  the  device  and  a  decrease  in  the  beam 
diameter.  In  the  binary  diffractive  element,  this  has  been  expressed  by  a  concomitant 
increase  in  the  frequency  of  the  chirped  grating.  Thus  for  best  results  the  diffractive 
element  should  be  designed  for  the  actual  beam  diameter  in  use.  This  would  also 
minimize  the  undesired  beam  fanning  inherent  in  the  use  of  a  chirped  grating,  as 
discussed  earlier.  The  perfonnance  of  the  diffractive  device  should  be  adequate  for  many 
applications  and  replicated  devices  could  be  fabricated  at  low  cost. 
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Appendix  A 

The  Taylor  series  expansion  of  the  index  of  refraction  n  of  the  refractive  element,  around 
the  point  of  symmetry  (x=0)  can  be  written  as 


/  \  un 

?(x  )=  n0  +  al  — x  +  a2  — 
dx  "  dx 


dn  d2n  ,  d3n  ,  d4n  4 


-x"  +a 


3  a  3  X  a4  a  4 


dx3 


dx4 


x  +• 


(Al) 


where  no  is  the  index  of  refraction  at  the  origin  (i.e.  the  point  of  symmetry),  and  ai,  ai ,  <73, 
<74  etc.  are  the  Taylor  expansion  coefficients  given  by 
_i  _  1  _  1  _  1  _  1 

CL  \  X  ■)  ^  ^  «.  CL  ^  «.  CL  a  »  CL  c  it  CL  /  • 

2!  3  3!  4!  5!  6! 

n  is  a  symmetric  function  around  the  origin.  This  causes  all  the  odd  powers  of  x  to  be 
cancelled  in  the  expansion  and  we  just  keep  the  even  powers. 


/  \  d^n  2  d4n  4 

n[x)  =  n0  +  a2 — 2x  +o4 — jx  -\ — 
dx  dx 

Similarly,  the  Taylor  series  expansion  around  the  coordinate  x+Ax  is 


(A2) 


n(x  +  Ax)  =  n0  +  a ,  (x  +  Ax)2  +  aA  (x  +  Ax)4  +  • 

dx  dx 


(A3) 


This  leads  to  the  variation  in  the  index  of  refraction  to  be, 
Aii(x)  =  ii(x  +  Ax)  -  n{x) ; 


(A4) 


since,  < 


(x  + Ax)2  -x2  =  A(x2)  =  2xAx 
(x  + Ax)4  -x4  =  A(x4)  =  4x3Ax 
(x  + Ax)6  -x6  =  A(x6)  =  6x5Ax 


we  have. 


dn 


An  =  2a-,x — 3- Ax +  4  aAx  . 

z  2  *+  n  4 

OX  ox 


d4n  .  ,  s  d6n 


Ax +  6 a,x  :L^-Ax, 

O  no  7 

ox 


(A5) 


which  is  Eqn.  (17)  in  the  body  of  the  paper. 
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We  consider  the  beam  retardation  for  a  beam  of  width  AD  traversing  a  graded  index 
material  whose  refractive  index  is  described  by  n(y).  The  value  of  A n  in  Eqn.  (A5)  can 
now  be  substituted  in  the  expression  for  beam  retardation  in  our  design  equation  (18),  so 
we  have: 

(A6) 

or,  after  the  substitution  and  dividing  both  sides  by  /,  we  have: 


-^-(An)(AD)l  =  — ^-sin-1 
A, 


t  X 

y 


m 


v„ 


2  n 

T 


^4  C\6 

on  A  3  o  n  ^  s  o  n 


lai  —< 

oy  oy 


■  +  6  a6y 


8y 


(AD)  =  ^sin  1 
col 


f  \ 

y_ 

vToy 


-A 


y 

Z2l 


Multiplying  both  sides  by  y  and  dividing  both  sides  of  the  equation  by  2ti(AD)/a 
we  obtain: 


6  86n  Agy  .  , 


2  d2n  4  d4n 

2a2y  ^+4«4v  tt+Kj  -  6  0  , 

oy  oy  oy  2  ncol 


sin 


f  \ 

y 


y  o) 


A 


y 


(AD)2tt  Z2l 


We  now  substitute  ao  =  2a2,  ai  =  4a4  and  a2  =  6a6  and  get: 


2„  a4„  a6„  a_.  f 


2  8  n  4  d  n  6  8  n  Agv  .  , 
ao y  +  +  — ^  =  — ^rSin 


a/ 


5/ 


d  y  2  tzcoI 


_y_ 

vWy 


y 


■— —  =  F(y) 
(AD)2tt  Z2l 


This  is  a  homogeneous  differential  equation  of  the  form 


ao  y 


dpn 


_!  dp-ln 


dvl 


-  +  ay 


dn 


- an  ,  v  —  +  an  =  F(y )  . 

dyp-1  dy  p 


(A7) 


(A8) 


(A9) 


(A10) 


One  way  to  solve  this  differential  equation  is  to  use  the  Cauchy-Euler  approach, 
through  variable  substitution,  where  we  substitute  y  =  e. 

Substituting  v  =  e  -A  t  =  In y,  for  y  >  0  in  the  equation  yields 


.  d6n  ,  din  ,  d4n  ,  d3n  ,  d2n  ,  dn 

An  7~  +  A,  h  An  —  A  4-,  h  vt4  r-  +  Ar  -  =  Cir(t ) 

dt6  dt 5  -  dt 4  3  dt2  dt 


(All) 
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where, 


G(t)  =  B0e‘ sin-1 - Bxe ", 

It<J 

and 

4  =  —,4  =  --,4  =  -—,4  = 

120  8  120  3 

B  _  AS^  B 
0  27t(AD)ojI  ’  1  2tt(AD)Z21' 


(A  12) 


49  29 

- ,4  =  — ,4  =  -3, and 

40  4  20 


This  equation  is  not  trivial  to  solve.  We  expanded  the  arcsine  function  in  a  series  to 
the  seventh  order: 


i  3 1  o  5 1  i  r  7 

1  e  3  e  15  e 


To )  To  6  To  40  To  336  To 


(A  13) 


This  yielded  the  expression. 


,  d6n  .  d5n  ,  d4n  ,  d4n  .  d2n  ,  dn 

4  — r  +  4  — r  +  4  — t  +  4  — r  +  4  — r  +  4  — 

0  dt 6  1  dt5  -  dt4  3  dt4  4  dt2  5  dt 

2 1  i  At  o  6t  i  r  „&t 

To  6  yd  40  y05  336  y07 


(A  14) 


Even  with  this  simplification  it  is  not  possible  to  find  an  exact  solution  to  the 
problem.  Therefore  we  considered  simplified  solutions  for  second  order  and  fourth 
order  expansions  of  the  equation  (A9). 

(a)  Second  order  solution 

Through  ignoring  the  4th  and  6th  order  differential  terms  in  Eqn.(A9),  by  following 
the  same  steps  as  we  did  in  obtaining  equation  (A  14)  which  included  these  terms  to 
the  6th  order  in  its  derivation,  we  found  the  following  relationship: 


d2n  dn 


an  „  ,  .  _i  e  „ 
- =  Bne  sin  —  -  B,e 


dr  dt 


(A  15) 
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The  solution  of  the  above  differential  equation  is: 


,  ,  5 By'  B/‘  B0e4'  e21  (Bq  -  Bx  v0 )  _ 

n(t)  = - 5 — t  +  — ^ — T  +  -5— r  +  — — - — +  C,  +e'C 

—  y*.  /  a  r\  r\  i  rn  r\  S  r\  V  L 


6272 V  400  v05  72  Vq 


2  To 


Substituting  v  =  e  yields: 


-  w = + ^4 + ^4 + y  2  ( V g|  y° } + Ci + ^ 


6272  v0  400v’()  72v0J 


2v0 


(A  16) 
(A  17) 


Our  assumption  since  the  beginning  was  that  n  should  be  symmetric.  Therefore  C2 
should  be  zero,  so  that  we  may  have  a  symmetric  function.  Through  multiplication  of 
both  the  numerator  and  the  denominator  of  the  terms  in  Eqn.  (A  17)  by  vo,  assuming 
C2  =  0,  and  rearranging,  it  is  possible  to  obtain: 


n{y) 


_  5B0y0 

6272 


_y_ 

y0) 


5oTo 


400 


z_ 

yoj 


B0y0 


72 


_y_ 

y0) 


+ 


(5o-AToK 


z_ 

y0J 


+c, 


1 9 


(Ai8) 


which  is  the  2nd  order  solution,  shown  in  equation  (19)  in  the  main  text. 


(b)  Fourth  order  solution 


Considering  terms  to  fourth  order  on  the  left  side  of  Eqn.  (A9),  the  equation  reduces 
to: 


2  8  n  4  d  n 
ao y  irr  +  W  —  j  =  F{y) 


8  y 


df 


(A19) 


Following  the  same  approach  as  before,  the  4th  order  equation  can  be  written  as: 


1  d4n  d3n  17  d2n  _  dn 

- - r  + - 5 — 2 —  =  By  sin 

6  dt4  dt 3  6  dt"  dt 


r e " 


vToy 


-By 


it 


(A20) 


which  has  the  following  solution: 
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=  5J°-V‘,+  By^+/(Ba-6B,yi,)+c  c 

225792  v07  7200  Vq  576v03  12v 


V23 

2 


In  v  +  C4sin 


V23 

2 


(A21) 


For  n  to  be  an  even  function, 

C2=  0,  C3  =  0,  C4  =  0 ,  so  we  have 

5Buy0  (  y  V  t  B0y0  f  y  Y  [  B0y0  f  y  V  |  y2 (B0  -  6Bxya )y0  ( +c  (A22) 

•  225792 ^ jo J  7200 {y0)  576  [yj  12  [yj  *’ 

which  corresponds  to  equation  (20)  in  the  main  text. 
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Figure  Captions: 

Fig.l.  Ray  Diagram  of  Resonant  Scanner  with  Refractive/Diffractive  Corrective  Element. 
Fig  2.  Schematic  diagram  illustrates  wavefront  bending  due  to  phase  lag  between  the  two 
sides  of  the  beam  as  the  beam  propagates  in  a  graded  index  of  refraction  material. 

Fig  3.  Plots  of  sin'x  and  its  polynomial  expansion  to  the  seventh  and  the  ninth  orders. 
The  solid  line  is  the  plot  for  sirf'(x),  the  dotted  line  the  polynomial  expansion  for  the 
seventh  order.  The  dashed  line  is  the  polynomial  expansion  to  the  ninth  order 
Fig.  4.  Linearized  Scanning  Geometries  with  FOV,  greater  than,  equal  to  or  less  than  that 
for  the  sinusoidal  scanner. 

Fig.  5.  Diffractive  element  for  the  second  order  solution  with  (a)  =  100,  vo=Tcm  and 

AD 

AD=0.01  cm,  b)  -^-  =  200,  vo=Tcm  AD=0.005  cm,  (c)  -^-  =  500,  vo=lcm  AD=0.002 
AD  AD 

cm,  as  generated  in  MATLAB.  The  calculated  grating  spacing  is  a  smooth  chirp.  The 
slight  banding  seen  in  the  figures  is  an  artifact  of  the  way  MATLAB  displays  these 
patterns. 

Fig.  6.  Plots  of  the  graded  index  of  refraction  as  a  function  of  the  normalized  distance 
from  the  optical  axis  for  the  second  order  solutions,  (a)  For  constant  vo  and  varying  beam 
diameters  AD  (b)  For  varying  vo  and  constant  beam  diameter  AD. 

Fig.  7.  Diffractive  element  for  the  fourth  order  solution  with  (a)  =  100,  vo=Tcm  and 

AD 

AD=0.01  cm,  b)  -^-  =  200,  vo=lcm  AD=0.005  cm,  (c)  -=^-  =  500,  vo=lcm  AD=0.002 
AD  AD 

cm,  as  generated  in  MATLAB.  The  calculated  grating  spacing  is  a  smooth  chirp.  The 
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slight  banding  seen  in  the  figures  is  an  artifact  of  the  way  MATLAB  displays  these 
patterns. 

Fig.  8.  Plots  of  the  graded  index  of  refraction  as  a  function  of  the  normalized  distance 
from  the  optical  axis  for  the  fourth  order  solutions  (a)  For  constant  v0  and  varying  beam 
diameters  AD  (b)  For  varying  yo  and  constant  beam  diameter  AD. 
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